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Abstract 

The stability/instability condition of a turbine rotor with axisymmetric supports is determined in 
the presence of gyroscopic loads and rub-induced destabilizing forces. A modal representation of 
the turbine engine is used, with one mode in each of the vertical and horizontal planes. The use 
of non-spinning rotor modes permits an explicit treatment of gyroscopic effects. The two 
linearized modal equations of motion of a rotor with axisymmetric supports are reduced to a 
single equation in a complex variable. The resulting eigenvalues yield explicit expressions at the 
stability boundary, for the whirl frequency as well as the required damping in the presence of the 
available rub-induced de-stabilization. Conversely, the allowable destabilization in the presence 
of the available damping is also given. 

Introduction 

Simulation of blade to case rubbing is critical to the complete understanding of the dynamic 
behavior of rotordynamic systems such as aircraft engines. Aircraft engines experiencing 
unbalance loadings, blade-outs, or mechanical misalignments may incur blade to case rubbing 
which may adversely affect the dynamic behavior of the engine and supporting structures. In 
some situations, the rubbing may become excessive and may even lead to structural instability 
and in severe cases catastrophic failure of the structure. 

The purpose of the present paper is to determine the conditions for stability/instability of a 
turbine engine with rub-induced forces in the presence of gyroscopic forces. To develop a direct 
determination of the stability/instability conditions, the equations of motion with rub forces 
obtained in Reference [1] are used as the starting point for the present investigation. In this study, 
the vertical and horizontal planes are each represented with one non-spinning vibration mode. 
Stability/instability is then investigated by examining the corresponding eigenvalues of the 
characteristic equation. 

The thrust of this work is with a rotor with axisymmetric support stiffness. This model allows for 
a complex-variable representation of the rotor system, and a more direct determination of 
conditions at the stability-instability boundary. The notion of neutral stability at the boundary 
permits obtaining two real equations from the complex eigenvalue problem. Equations are 
derived for the possible whirl frequencies and the required damping (or the allowable de- 
stabilization), at the stability/instability boundary. 

Blade Rub Forces and Moments 

Limited data and experience with mb tests on turbofan engines and test rigs indicate that blade 
tip rub forces, resulting from blade tip and case interactions, are in a direction that is 
perpendicular to the blade tip chordline. This phenomenon is similar to the behavior of an airfoil 


traveling through a fluid where the airfoil motion generates a lift and drag forces. From tests, the 
tangential mb blade force was found to have an axial component that would move the blade 
forward (or aft), and a circumferential component that opposes the rotor rotation. The 
hydrodynamic analogy expresses the blade mb force as proportional to the blade tip’s mean 
velocity and the instantaneous angle of attack or incidence. It differs from the typically used 
friction model of mb in that the friction model only generates a force in the direction opposing 
the rotor spin (i.e., drag) while the plowing model generates this drag force, as well as a “lifting” 
force that is perpendicular to the direction of rotor spin. A more detailed description of this mb 
model is provided in Reference [1]. 

An analytical expression, analogous to the dynamic pressure force from a change in momentum 
is used to develop the plowing model. This model is consistent with expressions of impact forces 
and hydrodynamic pressure. This was necessitated by the observed directionality of the blade tip 
“tangential” mb forces, i.e., perpendicular to the blade chordline. 

The hydrodynamic force is defined as: 

f =Dv 2 R 2 Yo A=DV 2 y 0 A (1) 


where f is the total tangential force on the blade from rubbing, V is the blade tip velocity, Yo * s 

the blade chordline angle of incidence in the rotating frame, and A is the frontal area of the 
robbing material (i.e., chord length x mb depth). D is a coefficient that must be determined from 
test or some other analysis procedure. It is most likely a function of several factors such as the 
material’s Young and Shear moduli, density, hardness, machinability, feed depth, dynamic shear 
strength, smoothness and probably others. 

The loads in the fixed coordinates system are: 


f x = f sin Yo sinjb - K r x 

(2a) 

f y = -f sin Yo cos|5 - K r y 

(2b) 

f a = -f cosYo cosf>R 

(3a) 

fg =f cosYo sin f>R 

(3b) 


Where Kr is the radial stiffness. Noting that: 

sinji = and cosf> = 


And substituting into (2) and (3), the rub loads become, 


f x =D4/ 2 R 2 Y 0 ysinY 0 -K r x 

(4a) 

f y = -D<j/ 2 R 2 Y 0 xsinYo -K r y 

(4b) 

f a =-Dij/ 2 R 3 Y 0 xcosYo 

(5a) 

f e =Dq/ 2 R 3 Y 0 ycosY 0 

(5b) 


Modal TVansformation of the Rub Loads 

As a first step to simplifying the equations of motion, the translational and rotational rub forces 
in (4) and (5) are transformed to modal coordinates. The following considerations are used to 
obtain the modal rub forces: 

1. There is one modal function in each of the horizontal {x, a) and vertical (y, 0) planes. 

2. The rotational (a, 0) deflections are the slopes of the bending displacements. 

3. Because of the {x, y, a, 0) sign conventions, it is seen that: 

- dy , dx. 

dz dz 


4. The modal functions in the fixed axes are non-rotating. 

5. The modal functions are mass normalized. 

Now let the modal displacements be as follows: 

x = xh(t) (6a) 

y = yv(t) (6b) 

a = ah(t) (7a) 

0 = -0v(t) 


(7b) 


Where x , y , a , 0 are mode shapes and h (t) and v(t) are horizontal and vertical plane 
displacement generalized coordinates, respectively. 

Following the Principle of Virtual Work, the modal forces in the horizontal and vertical planes 
are: 

h = J(xf x +af a )ds (8a) 

f v = {(yf y -0fe)is (8b) 

Where the rub forces are integrated over the tangential extent of the rub, “s”, then summed over 
all rubbing blades. 

Making the physical-to-modal displacement substitutions, the modal rub forces become: 

f h = j{Dtj> 2 R 2 Yo(xysinYov(t)-RaxcosYoh(t))-KrX 2 h(t)}ls (9a) 

f v = J {- Dij/ 2 R 2 Yo (xy sin Y 0 h(t) + R0 x cos Y 0 v (t))~ K r y 2 v(t)jds (9b) 

Defining the following modal integrals: 

K x = |K r x 2 ds 

S xy = fxy ds 
S 0y = J©yds 

K y = jK r y 2 ds 
S ax = Jaxds 

Then substituting into (9a) and 9(b) yields: 

fh = Dtp 2 R 2 Y 0 {s xy sin Y 0 v(t)-RS ax cosYoh(t)}— K x h(t) 


(10a) 


f v = Dv 2 R 2 Y 0 {-S xy smY 0 h(t)-RS ex cosY 0 v(t)}-K y v(t) 


(10b) 


The following variables are defined to ease the required algebraic manipulations: 

S 0 = D(j/ 2 R 2 Y 0 S xy sin y 0 

S h =Dv 2 R 3 Y 0 S ay cosYo 

S v =Dv 2 R 3 Y 0 S ey cosYo 

Substituting into (10a) and (10b) the modal mb force equations are: 


fh =S„v(t) — K x h(t)-S„h(t) 

(11a) 

fv =-S 0 h(t)-K,v(t)-S,v(t) 

(lib) 


The Modal Equations of Motion 

The modal equations of motion are obtained by combining the modal mb forces, (11a) and (lib), 
with the modal inertia, stiffness, damping and gyroscopic forces. The latter forces are obtained 
by using the translational and rotational values of the horizontal and vertical modal vectors to 
transform the physical inertia, stiffness, damping and gyroscopic properties to modal 
coordinates. Performing these transformations and using the modal mb forces derived in the 
previous section yields the following modal equations of motion (see also References [2] and 

[3]): 


ii + u^h + 2u; h £ h h + 2q»gv = f h (12a) 

v + o» 2 v + 2 uj v £ v v - 2q»gh = f v (12b) 

Substituting the rub forces and collecting coefficients of the same dependent variables we have 
the following: 

ti + (u/h +K X +S h jh + 2ui h 5 h h + 24>gv-S 0 v = 0 (13a) 

v + (u ; 2 + K y + S v ]v + 2ut v £ v v -2qjgh + S 0 v = 0 


(13b) 


In general, the system is non-axisymmetric in the rotor/support modal stiffness, as well as in the 
damping. Also, unless the horizontal and vertical modal displacements are equal, the direct mb 

force coefficients Sh, and S v , K x and Ky , also are not equal. However, many incidents of turbine 
engine whirling exhibit circular or almost circular orbits. Furthermore, one may also be able to 
construct highly elliptical orbits from identical uncoupled horizontal and vertical mode shapes 
even when the corresponding modal stiffnesses or frequencies are grossly different. Therefore, to 
simplify the modal equations of motion it will be assumed that the horizontal and vertical mode 
shapes are the same. However, the support stiffnesses and modal frequencies may be unequal. 

With this consideration the following additional simplification is made: 

s h = s v = Sj 

K x = K y = K 


Where the local case and rotor radial stiffness is assumed to be isotropic and axisymmetric. 
Additionally, the modal damping is assumed to be identical in both planes. The damping 
symmetry is justified since damping values are difficult to obtain and in most analyses is an 
educated guess, at best. Assuming symmetric damping: 


l = 2Wh€h = 2ut v £ v 


Finally, the modal equations of motion are: 

h + + K + Sj |h + £h + 2ijigv — SqV = 0 (14a) 

v + (mJ + K + S 1 )v + Cv-2vgh + S 0 h = 0 (14b) 

Where the support non-axisymmetry is reflected in the horizontal and vertical modal frequencies, 
ut h and w v . 

Stability of Rotor with Axisymmetric Supports 

It had been shown by various investigators such as Gallardo [4], Smith [5], Ehrich [6], Alford [7] 
and others, that the least stable rotating system is axisymmetric. Thus one may initially perform a 
stability analysis of a rotor with an axisymmetric support. Then, if the axisymmetric system is 
stable, no additional analysis may be required, since a non-axisymmetric supported system will 
always be more stable than the identical rotor with symmetric supports. 

In the special problem of an axisymmetric rotor system the stability/instability conditions are 
obtained with greater simplicity and directness. The modal equations of motion of rotor with 
axisymmetric support stiffness are obtained from the previous section by making the modal 
frequencies equal. Thus: 


(15a) 


h + (u>o +K + S!)h+Ch + 2ij)gy-Sov = 0 
v + (u/g + K + S 1 ]v + £v-2tt>gh+S 0 h==0 (15b) 


Where: w 0 =w h =w v . 


Equations (15a) and (15b) are combined into a single equation by introducing the complex 
variable z, where z = h + iv . Multiplying the second equation of motion, (15b), by (i) and adding 
it to the first, (15a), and using the complex variable definition, we have 

z + (u)q +K + S 1 )z + ^z-i2^gz + iS 0 z = O (16) 

Assuming a solution in the complex exponential form: 

z = z 0 e iAt = z 0 e i(A+iB > = z 0 e~ Bt e iAt 

Where the coefficient Zo and the eigenvalue A are complex quantities and “A” and “B” are real 
numbers. 

Furthermore: 

Real (A ) = A = frequency (real quantity) 

Imag ( A ) = B = attenuation or amplification factor (real quantity) and: 

B > 0, stable 
B < 0, unstable 

Stability/Instability Boundary of the Axisymmetric System 

Substituting the exponential solution into the complex variable equation of motion (16), one 
obtains: 



-2Ati>g-(tUg +K+S 1 ^-i(A£ + S 0 )jz = 0 


(17) 


The condition for the existence of a nontrivial solution requires that the polynomial coefficient of 
“z” must vanish: 


Ft*) = |* 2 - 2Al$ig - (wi + K + S, | - i(A? + S„ )}= 0 


(18) 


At the stability-instability boundary, the solution is simple harmonic, which implies that A = 
Real at the boundary. Therefore, since all the parameters in the Eigen polynomial, (18), are real 
and non-zero, then the following holds at the stability/instability boundary: 

Real part of F (A) = 0 

Imaginary part of F ( A ) = 0 

And we have two equations at the stability boundary: 

A 2 -2Aii»g-(uu 0 2 +K + S 1 ) = 0 (19a) 

A£ + S 0 = 0 (19b) 

Equation (19a) is solved for A which gives the possible frequencies at the boundary. 

A = (pg±V<i> 2 g 2 + (uto +K + SJ 

Using the frequency, A determined above and Equation (19b), the damping required for stability 
is: 

So 

A 

Since the term So is always destabilizing (So> 0), for a system to be stable, the system damping 
£ must be positive (£ >0). Therefore: A = A. < 0, at the boundary because B = 0 at the 
boundary, and £ and Sq are positive. Furthermore, A = A + > 0, is always stable because if A > 0 
then B > 0. 

From the above, the results for the rotor-rub stability /instability boundary are: 

Whirl frequency: A_ = yg - -Jip 2 g 2 + (iUo +K + Sj) 



And the required damping, = - 


So 

A_ 


Comparison of the available damping to the required damping (in the presence of a de- 
stabilization) gives a measure of the stability of the system: 


C available ) £req 5 stable 
C available ( Creq > unstable 


Recall that So and Si contain the plowing rub coefficients, comprised of the rub constant, rotor 
speed, blade tip angle and tip radius and k r is the radial-rub restoring or snubbing force, which 
may or may not be dependent on the rub parameters. Recognizing: 


cot Y 0 = 


cos Vq 
sinYo 


Si may be written in terms of So as: 


Sj = S 0 Rp, where p = 



cot Yo 


Making the above substitutions into the expressions for the whirl frequency and required 
damping, we have, at the stability instability boundary, the following results. 

Whirl frequency: A_ = ipg - ^/q/ 2 g 2 + (cUq + K + S 0 Rp) (20a) 

§ 

Required damping: = ^ Q (20b) 

Vg - VV g + (utf + K + S 0 RpJ 

The rotor shown in Figure 1 is used to demonstrate the use of Equation (20) for determining the 
whirl frequency and required damping for a known quantity of rub. The rotor consists of a 
massless cantilever shaft with a disk attached at the free end. The properties of the rotor are 
shown in the figure. The non-rotating frequencies and modes shapes are computed from a finite 
element analysis. The first modal frequency and mode shape is: 

w 0 = u/ h = uj v =3.071xl0 2 rad/sec 


(l \J’\. 1X0 \ 

\ a / \0 / \ 7.144 xlO -2 / 


Disk 


M = 100. Ibm 

l p = 20,000. Ibm-inch 2 

Gravity = 386. inch/sec 2 

Radius = 5 inch 

Blade Angle = 30 degrees v 

MS 


Rotation 

Vector 


Rotor Shaft 
E = 30.X10 6 lb/inch 2 
I = 20. inch 4 



f - — Rub Spring 

K = 3.x103 lb/inch 
Rub Length = 5 inch 


Figure 1. Cantilever Rotor 


Substituting the above values, and the values for blade angle, disk radius, rub stiffness and 
length, into Equation (20a) and (20b), the whirl frequency and required damping are: 

A = 0.2644 ip - /+ -y/.0699 <(/ 2 + 14.45 x 10 4 + 37.04Dqt 2 

109.6 Dty 2 
A 

Figure 2 shows the result of plotting the whirl frequency and required damping as a function of 
rotor speed for a variety of dynamic rub coefficients. The whirl frequencies associated with both 
the forward, A + and the backward, A_ mode are shown in Figure 2a, however only the backward 
whirl frequency is significant since the backward whirl frequency is used to calculate the 
required damping for a stable system. The required damping as a function of rotor speed for 
various levels of rub is shown in Figure 2b. As expected and depicted in the figure, the required 
damping increases with both rotor speed and mb force. 





Rotor Speed (Hz) 


Figure 2a. Whirl Frequency at Stability/Instability Boundary 




Rotor Speed CHzj 


Figure 2b. Required Damping at Stability/Instability Boundary 


Allowable Destabilization for an Available Damping 

The damping relation may be solved for the allowable destabilization So, in terms of the available 
damping. In this context, Equation (20b) is solved (with the available damping assumed to be 
known), for the destabilization that can be allowed that will make the initially stable system 
approach neutral stability or the stability/instability boundary. 

Thus the corresponding relation is: 


So = -A? 


And the converse expression for the allowable destabilization in the presence of the available 
damping becomes, 


^0 Allowable £ Available S + (ui 0 + K + So^i owable Rp) j 


( 21 ) 



Equation (21) is complicated since So is on both sides of the equal (=) sign. However, with some 
algebraic manipulation, the allowable destabilization can be obtained explicitly. Rearranging 
Equation (21) yields: 


S 


^Allowable 


gftRp-2q>g) 

2 



( 22 ) 


Noting that So occurs skew-symmetrically in the modal equations of motion [Equation (13)], and 
that it has a positive sign (+) in the vertical modal equation and a negative (-) sign in the 
horizontal. Thus, in the mathematical context, only a positive (+) for So is admissible. The signs 
of So in the equations of motion reflect the unstable whirl direction, such that a positive sign in 
the vertical DOF (negative in the horizontal DOF) may indicate a backward whirl, while 
opposite signs would denote a forward whirl. This effect is also seen in the expression for the 
required damping. 

The assessment of system stability or instability can then be made by comparing the available 
destabilization to its allowable value, in the presence of the available damping. Thus: 


^Available ^ ^Allowable' s l a ^ e 


^Available ^0 Allowable' unstaWe 


Concluding Remarks 

The linearized modal equations of motion of a rotor with non-axisymmetric stiffnesses with both 
gyroscopic and rub forces have been obtained. The equations of the special case of a rotor with 
axisymmetric support have also been derived. The latter equations are used to determine the 
influence of the rub forces on rotor stability/instability boundaries. 

The two-mode representation (one in each vertical and horizontal plane) of the equations of 
motion was considered. For the axisymmetric support case, the modal equations of motion were 
reduced to one equation using a complex variable to describe the rotor motions. A method was 
derived to investigate the conditions at the stability/instability boundary. 

At the boundary, the whirl frequencies and required damping for the available value of the de- 
stabilizing coefficient are obtained. The converse, which is the allowable de -stabilizing 
coefficients in terms of the available damping, is also derived. These expressions are in closed 
form, and provide a convenient way to assess the stability/instability implications of both rub 
force parameters and the rotor system’s modal inertial, damping, elastic and gyroscopic 
properties. 

The analysis for the axisymmetric rotor system may be considered as a negative 
stability/instability criterion, in that an axisymmetric system is the least stable configuration. 


Thus, if this analysis indicates a stable rotor, an analysis of a rotor with non-axisymmetric 

support stiffnesses may not be required, since the latter configuration is always more stable than 

the former. 
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